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This article presents the compatibility of experimental data from neutrino oscillation experiments 
with a high-Am 2 two-neutrino oscillation hypothesis. Data is provided by the Bugey, Karlsruhe 
Rutherford Medium Energy Neutrino Experiment 2 (KARMEN2), Los Alamos Liquid Scintillator 
Neutrino Detector (LSND), and MiniBooNE experiments. The LSND, KARMEN2, and MiniBooNE 
results are 25.36% compatible within a two-neutrino oscillation hypothesis. However, the point of 
maximal compatibility is found in a region that is excluded by the Bugey data. A joint analysis of 
all four experiments, performed in the sin 2 28 vs Am 2 region common to all data, finds a maximal 
compatibility of 3.94%. This result does not account for additions to the neutrino oscillation model 
from sources such as CP violation or sterile neutrinos. 



INTRODUCTION 

Neutrino oscillations have been reported at three dif- 
ferent Am 2 scales: solar, atmospheric, and high- Am 2 . 
The solar [? ] and atmospheric [? ] best-fit results 
have been observed by several independent experiments, 
using various neutrino sources and techniques. The high- 
Am 2 result, from the LSND detector [? ], has yet to be 
reproduced. 

The LSND experiment observed a significant excess 
of events which are best fit by — > v e oscillations at 
the Am 2 ~ 1 eV 2 scale. The solar and atmospheric 
best fit results are at Am 2 ~ S^qa x !0~ 5 eV 2 and 
Am 2 ~ 2.4±^J x 1(T 3 eV 2 , respectively This wide 



spread of Am 2 scales cannot be accommodated by the 
three neutrino mass states of the standard model. The 
LSND result is uniquely incompatible with other oscil- 
lation observations, and if verified would demand exten- 
sions to the standard model in the neutrino sector [? ] [? 

]■ 

Prior to 2007, two experiments, KARMEN2 and 
Bugey, performed searches for oscillations in the region 
of oscillation parameter space probed by LSND. KAR- 
MEN2 [? ] conducted an accelerator-based —f D e 
appearance search. The Bugey [? ] reactor experi- 
ment probed for oscillations using v e disappearance. Nei- 
ther experiment found evidence for neutrino oscillations. 
However, a joint analysis [? ] between LSND and KAR- 
MEN2 found both were compatible with a two-neutrino 
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oscillation hypothesis at 64% confidence level (CL). for 
Am 2 > 0.2 eV 2 , in a region not covered by Bugcy. 

The MiniBooNE experiment, located at Fermi Na- 
tional Accelerator Laboratory, was designed to fully ex- 
plore the LSND result. In 2007, MiniBooNE published 
results from a — > v e appearance oscillation search [? 
] . MiniBooNE observed no significant excess of events in 
an energy range from 475 MeV to 3 GeV. MiniBooNE 
is presently collecting anti-neutrino data, for use in an 

— > D e appearance oscillation search. 

This analysis presents results from the combination of 
LSND, MiniBooNE, KARMEN2, and Bugey. It is moti- 
vated by a need to determine if the LSND excess may be 
the result of two-neutrino oscillations, in light of these 
three null experiments. Results presented in this arti- 
cle make use of the MiniBooNE neutrino data set, and 
do not include the unpublished anti-neutrino data. The 
compatibility found in this analysis is valid within the 
framework of standard 2-neutrino oscillations. 
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Using this pcr-cvent probability in a likelihood fit, the 
log likelihood becomes 



(xj - f(xi\a)) 2 
2a? 



^]lnV2^cr 4 . (3) 



The second sum in the ln£ equation does not typically 
depend on the theory parameters a. Equation [3] is mini- 
mized by minimizing the familiar x 2 function, 



(xi - f(xi\a)f 



(4) 



INPUT DATA 

Data from each experiment are provided in a two- 
dimensional (2-D) grid of sin 2 29 vs Am 2 . The value 
at each grid point represents the agreement between the 
observed data and a two-neutrino oscillation hypothesis, 
with a signal appropriate to the oscillation parameters at 
that point. The data sets come in several different for- 
mats (log likelihood (ln(L)), Aln(L), x 2 )i spanning dif- 
ferent sin 2 29 vs Am 2 ranges. An optimal compatibility 
calculation would make use of the absolute x 2 > as °P~ 
posed to the Ax 2 , which is the change in \ 2 between each 
grid point and the experiment's best fit point. However, 
we were unable to obtain the absolute % 2 information 
from all input experiments. Therefore, our compatibility 
calculation can only make use of relative A\ 2 informa- 
tion. All input data is transformed into a A% 2 surface 
in sin 2 29 vs Am 2 space, with common sin 2 29 vs Am 2 
binning. 

The transformation from the input ln(L), Aln(L) grid 
to a A% 2 grid is derived as follows. In the technique 
of maximum likelihood fitting, a per-event probability 
p(xi\a) is constructed, where the event-measured 

quantities and a are the theoretical parameters. The 
goal is to maximize the likelihood, C(a), the probability 
of all events in the sample, assuming the given model 
probability p(xi\a) for each event, 

C{a)=n iP {xi\a). (1) 

The technique of x 2 fitting is a special case of likelihood 
fitting in which the per event probability is a Gaussian 
distribution, 



From the point of view of minimization, contours, and 
interpretation of results, there is an equivalence between 
the likelihood and the x 2 functions, given by 

-2lnC(a) = x 2 (a). (5) 

To convert between the input ln(L) or Aln(L) data and 
the x 2 or Ax 2 use d in this analysis, the input data is 
multiplied by a factor of -2. The validity of this conver- 
sion technique has been verified by comparing calculated 
allowed regions found using the Aln(L) and A% 2 grids, 
for LSND and KARMEN2, with those published by these 
experiments. 

The input experiments published observation and limit 
curves using two different methods. The first method is 
a two-dimensional (2D) global scan that calculates the 
Ax 2 or Aln(L) with respect to the global best fit point 
across the entire grid. LSND and KARMEN2 calculated 
their results using this method. The second method is 
a one-dimensional (ID) raster scan that calculates the 
change with respect to the local best fit point in each 
Am 2 row. The Bugey and MiniBooNE experiments used 
this method to produce their exclusion curves. Given the 
mixture of methods used to report results from the input 
data, we have performed our compatibility calculation 
using both methods. For each input experiment we create 
two A% 2 grids - one using the 2D global scan method, 
and one using the ID raster scan method. 

The LSND data are provided as a 2D histogram of 
ln(L) values containing the decay-in-flight and the decay- 
at-rest results. The input grid covers 0.000313 to 1.01 in 
sin 2 29 and 0.0098 to 101.16 eV 2 in Am 2 . The LSND 
ln(L) grid is first converted into a Aln(L) grid, and then 
multiplied by -2 to produce a A% 2 grid. The conversions 
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from ln(L) to Aln(L) and Aln(L) to A% 2 are tested by 
calculating the 2D 90% and 99% CL allowed regions by 
stepping away from the global best fit point a certain 
number of units in Aln(L), A\ 2 space. Aln(L) units are 
2.3 (90% CL) and 4.6 (99% CL); A X 2 units are, for a 2D, 
2 degree of freedom scan: 4.61 (90% CL) and 9.21 (99% 
CL). Both tests properly reproduce the published LSND 
result [? ]. 

The data from the KARMEN2 experiment are Aln(L) 
values covering a range of 0.000316 to 1 (sin 2 29) and 

0. 01 to 100 eV 2 (Am 2 ). Each point is multiplied by -2 
to produce a A\ 2 grid. A cross-check is performed us- 
ing the A% 2 grid, where the probability is calculated at 
each point on the grid using 2 degrees of freedom (DOF) . 
Points where the probability crosses 10% delineate the 
90% exclusion band. This test correctly finds the 90% 
CL exclusion band from the KARMEN2 publication [? 

]■ 

We were unable to obtain data directly from the 
Bugey collaboration. However, a recent global analysis 
of Gallium and reactor v e disappearance data describes a 
method used to reproduce the x surface of Bugey, com- 
plete with full systematic errors [? ] . The authors kindly 
provided us with their full x 2 surface (sin 2 29: 0.01 to 

1, Am 2 : 0.01 to 100 cV 2 ). A cross-check was also per- 
formed on this data by applying the raster scan method 
and stepping away from each local best fit point by 2.71 
Ax 2 units to find the 90% CL exclusion band. Using this 
method we are able to reproduce the published Bugey re- 
sult [? ]. 

The MiniBooNE data is expressed as a x 2 format from 
0.0001 to 0.4108 in sin 2 26 and 0.0488 to 51.13 cV 2 in 
Am 2 . This analysis utilizes the MiniBooNE data from 
475 MeV to 3 GeV in neutrino energy; the low energy 
region (below 475 MeV) is not considered in the compat- 
ibility calculation. 

COMPATIBILITY CALCULATION 

The compatibility calculation uses a method developed 
by Maltoni and Schwetz [? ] to answer the specific 
question, ''How probable is it that all experimental results 
come from the same underlying two-neutrino oscillation 
hypothesis?". First, a A\ 2 grid is constructed for each 
experiment as described in the previous section. The 
individual grids arc then summed together to produce 
one summed Ax 2 grid. The compatibility test statistic, 
Xmini i s the minimum of the summed Ax 2 grid. Xmin 
lows a x 2 distribution with P c degrees of freedom, where 
P c is the sum of the total number of independent pa- 
rameters minus the number of independent parameters 
estimated from the data. For example, the combination 
of the 2D MiniBooNE and 2D LSND results yields four 
total independent parameters; each experiment indepen- 
dently measures sin 2 20 and Am 2 . Two parameters are 



estimated from the data (sin 2 29, Am 2 ), resulting in a P c 
of 2. The final compatibility is the x 2 probability of Xmin 
using P c degrees of freedom. In the analysis using two 
experiments P c is 2, for three experiments P c is 4, and 
for all four experiments P c is 6. 

This method is designed to be robust against cases 
where the x 2 minima of the individual data sets are very 
low, and when several parameters are fitted to a large 
number of data points. It reduces the problem that a 
possible disagreement between data sets becomes diluted 
by data points which are insensitive to the crucial pa- 
rameters. 

Of course, there are limitations to this method. This 
method does not take into account the absolute good- 
ness of fit of each individual experiment at its own best 
fit point. It is also valid only for truly statistically in- 
dependent data sets. Theoretical uncertainties in similar 
experiments may introduce correlations between the vari- 
ous results; for example, LSND and KARMEN2 have the 
same neutrino beam energy spectrum which may result 
in similar neutrino interaction errors. However, a pre- 
viously reported combined analysis of LSND and KAR- 
MEN claims the two experiments may be considered in- 
dependent [? ] . MiniBooNE and Bugey are not expected 
to have any uncertainties in common with the other ex- 
periments. 

ALLOWED REGION CALCULATION 

Combinations of experiments which result in a com- 
patibility of greater than 10% are further explored to 
locate any remaining allowed regions. Allowed regions 
are indicated by closed contours in the sin 2 29 vs Am 2 
plane. Contours which do not close form exclusion bands; 
parameter values situated to the right of the bands are 
excluded at a given CL (typically 90 and 99% CL). The 
allowed regions indicate where the oscillation parameters 
would lie, at a given confidence level, assuming all exper- 
imental results can arise in a framework of two-neutrino 
oscillations. The calculated compatibility is the metric 
for how valid this assumption is. 

The allowed regions calculation follows the prescription 
of Roe [? ]. Each experiment's Ax 2 grid is converted into 
a Ax 2 probability grid, using an appropriate number of 
DOF (two for the global scan analysis, one for the raster 
scan analysis) . The final combined probability at a given 
point can be obtained from the product of the individual 
probabilities, x. The result is a sum of powers of the 
absolute value of the logarithms of x, 

n—l -j 

Prob = x- y"--\ln j (x)\ (6) 
j=o J ■ 

, where x is the product of the individual probabili- 
ties and n is the number of experiments being included. 
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Points where the probability crosses 10% bound the 90% 
confidence level allowed region; points where the proba- 
bility crosses 1% bound the 99% region. 

There has been much discussion regarding the number 
of DOF that one can use with the \ 2 grids in the allowed 
region calculation, and how the DOF changes across the 
grids [? ] . We examined the change in DOF across the 2D 
grid, using the Feldman-Cousins frequentist method [? ]. 
Our study of the DOF finds that, in general, use of 2 
DOF is valid across the 2D grid. This breaks down for 
points with high Am 2 and high sin 2 28 values: Am 2 > 
10 eV 2 and sin 29 > 0.01. However, the approximation 
of 2 DOF is still valid for sin 2 20 > 0.01 and Am 2 < 10 
eV 2 . The region where the 2 DOF approximation is no 
longer valid is an area which does not contain any allowed 
regions from LSND, and as such should not impact this 
analysis. 

The compatibility is defined as the A% 2 probability at 
the best fit point. The combination of experiments re- 
duces the number of independent parameters used in the 
probability calculation. The Maltoni-Schwctz method 
can easily accommodate the change in degrees of free- 
dom [? ]. It is not clear how to include information 
about the number of degrees of freedom into the more 
traditional (Roe) method (Equation[6|). The Roc method 
applied as-is results in too high a compatibility, due to 
the inability of the method to consider a reduced degree 
of freedom. 

We have chosen to use the more traditional Roc 
method to calculate the allowed regions, and the Maltoni- 
Schwetz method to find the compatibility. Equivalcntly, 
we could have chosen to use the generic Maltoni-Schwetz 
method to find the allowed regions; both methods return 
identical values when evaluating the joint probability dis- 
tribution function. 

The range of sin 29 vs Am 2 common to all experi- 
ments is used for the compatibility and allowed region 
calculations. The Am 2 is restricted to 0.0488 to 51.13 
eV for all results. The sin 2 29 range, for results without 
Bugey, is 0.000317 to 0.4108. Results containing Bugey 
employ a sin 2 28 range of 0.01 to 0.4108. 

HOW TO INTERPRET THE RESULTS 

As previously discussed, the A% 2 grids used in this 
analysis are created in two ways. The 2D grids use the 
global best fit point internal to each experiment to pro- 
duce the Ax 2 . The compatibility extracted from the 
combined 2D A% 2 grid represents how probable it is that 
one could observe the experimental results if nature truly 
has two neutrino oscillations in this high-Am 2 region 
(0.0488 to 51.13 eV 2 ). This method also finds the most 
probable point for the true oscillations to exist across 
the evaluated phase space. The ID results use a raster 
scan method to find the local best fit point at each Am 2 , 



internal to each experiment, to produce the A^ 2 - The 
results from this method represent the compatibility at 
each Am 2 , if nature truly had two neutrino oscillations 
located at that Am 2 . 



RESULTS 

The 2D analysis reports a single value for the max- 
imum compatibility of the experimental data with the 
two neutrino oscillation hypothesis. These results are 
presented in Table [U The compatibility for the ID anal- 
ysis is a function of Am 2 , and is presented in graphical 
form in the following sections. All results are calculated 
with respect to an oscillation hypothesis valid in the Am 2 
region of 0.0488 to 51.13 eV 2 . This region has been fur- 
ther divided into three components ; low Am 2 indicates 
the region from 0.0488 to -1 cV 2 , medium Am 2 spans 
— 1 to ~7 cV 2 , and high Am 2 is >~7 cV 2 . These divi- 
sions are used to characterize the results in the following 
sections. 

We first present results from the combination of the 
three accelerator-based appearance oscillation experi- 
ments (LSND, KARMEN2, MiniBooNE). This combi- 
nation finds a high compatibility (25.36%) at low Am 2 
(Figurc[T|). The inclusion of the Bugey reactor disappear- 
ance data highly constrains the low Am 2 region, reducing 
the compatibility to a low level (3.94%, Figured]). In the 
second section, the compatibility and allowed regions are 
explored in various combinations of the null experiments 
(KARMEN2, MiniBooNE, and Bugey). In the third sec- 
tion, the result of 2.14% compatibility omits the KAR- 
MEN2 data. The final section discusses the combination 
of LSND and KARMEN2. 
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X X 32.21 1 0.066 1 0.4 



TABLE I: Maximum compatibility for a variety of combina- 
tions of the input experiments, found using the 2D A\ 2 grids. 
The last two columns indicate the Am 2 vs sin 2 26 location of 
the point of maximum compatibility. The X's indicate which 
experiments were included in the analysis. 



5 



LSND, KARMEN2, MiniBooNE 

First, we consider only results from the v e appearance 
searches. Figure Q] (top) displays the 2D A% 2 grid from 
the combination of LSND, KARMEN2, and MiniBooNE. 
The point of maximal compatibility (25.36%) is indicated 
by the star. The point of highest compatibility is not lim- 
ited by the sin 2 29 vs Am 2 grid boundaries, but is found 
in a region excluded by the Bugey data. The allowed re- 
gions for two-neutrino oscillations are shown in the bot- 
tom of Figure [TJ There are 99% allowed regions at low, 
medium, and high Am 2 . The only 90% allowed region 
is located at low Am 2 , which overlaps slightly with the 
LSND 90% allowed region. 

Compatibility values are also reported for the ID anal- 
ysis as a function of Am 2 . The maximum compatibility 
for the ID LSND, KARMEN2, MiniBooNE analysis is 
shown in the top of Figure [2] The ID analysis also finds 
a high compatibility (here almost 50%) at low Am 2 . In 
addition, this raster scan method allows for a ~25% com- 
patible region at medium Am 2 . Figure [2] (bottom) illus- 
trates the 90 and 99% allowed regions. While these re- 
gions appear to be shifted from the LSND signal region, it 
must be remembered that the LSND signal region shown 
is that found using a 2D analysis, not a ID scan. 

Finally, Figure [3] compares the LSND allowed regions 
found using the ID A% 2 method to those published by 
the LSND collaboration, found using the 2D scan. The 
ID allowed regions are quite large in comparison to the 
2D regions, and include areas which had no 90% allowed 
islands in the 2D scan. 
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LSND, KARMEN2, MiniBooNE, Bugey 

The inclusion of the Bugey data significantly changes 
the compatibility results. The Bugey sin 2 29 range has 
a lower bound at 0.01; analyses including Bugey are re- 
stricted to sin 2 29 of 0.01 to 0.4108. The combination 
of all four experiments has a 2D compatibility of 3.94%. 
Figure [4] shows the final 2D Ax 2 grid, for all four exper- 
iments. 

The ID analysis of all four experiments agrees quite 
well with the 2D results; the point of highest compati- 
bility is found in the low Am 2 region, and all results are 
no more than 5.2% compatible with having resulted from 
two- neutrino oscillations (Figure [5|). 

KARMEN2, MiniBooNE 

It is instructive to calculate compatibility and remain- 
ing allowed regions, in the absence of a positive LSND sig- 
nal. The 2D analysis finds KARMEN2 and MiniBooNE 
are 73.44% compatible with a two-neutrino hypothesis; 
there is a 73.44% chance that we would find these two 



FIG. 1: Top: Summed 2D Ax 2 compatibility grid from 
LSND, KARMEN2, and MiniBooNE. The star indicates the 
point of maximal compatibility (25.36%). Bottom: Allowed 
regions (90%, 99%) found for the 2D LSND, KARMEN2, and 
MiniBooNE joint analysis. Triangle points contain the 99% 
CL region, circle points contain the 90% CL region. The 
solid brown area is the LSND 90% allowed region; the solid 
light blue area is the LSND 99% allowed region. The vertical 
straight edge on the left arises from a sharp discontinuity in 
the LSND input grid. 

null results, in the presence of two-neutrino oscillations at 
these sin 2 26>, Am 2 values (0.147, 0.052 eV 2 ). However, 
the minimum in the A% 2 is outside of the region where 
either experiment has much sensitivity. For example, the 
compatibility at the lowest grid point (0.0003, 0.05 eV 2 ) 
still remains high at 53.52%. The summed A% 2 and al- 
lowed regions for oscillations are shown in Figure [SJ 

The ID analysis of KARMEN2 and MiniBooNE finds 
high compatibility that reaches 100% at low and medium 
Am 2 (top, Figure [7]). However, the 90 and 99% CL ex- 
clusion curves (bottom, Figure [7]) arc almost identical to 
those found in the 2D analysis. 
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Maximum 1-D Compatibility (LSND, KAR2, MB) 
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FIG. 2: Top: Maximum compatibility as a function of Am 2 
for the ID LSND, KARMEN2, MiniBooNE analysis. Bottom: 
Allowed regions for the ID LSND, KARMEN2, MiniBooNE 
analysis. 90% allowed regions exist in the low and mid Am 2 
regions. Triangle points contain the 99% CL region, circle 
points contain the 90% CL region. The vertical straight edge 
on the left arises from a sharp discontinuity in the LSND input 
grid. 



KARMEN2, MiniBooNE, Bugey 

If we ignore the positive LSND result, but now include 
Bugey, it is 27.37% probable that we would have found all 
three null results in a world with two-neutrino (Am 2 > 
0.0488 eV 2 ) oscillations. Please note that the point of 
maximal compatibility is limited by the boundary of the 
analyzed region. Figure [8] (top) shows the 2D compati- 
bility of all three null results. Figure [H] (bottom) presents 
the remaining allowed regions. The straight line on the 
left hand side is an artifact of the requirement that the 
analysis be performed over regions of phase space com- 
mon to all experiments. 

The ID analysis of KARMEN2, MiniBooNE, and 
Bugey (top, Figure [5]) produces a higher degree of com- 
patibility than the 2D analysis shown in Table [H but 
agrees with the remaining allowed regions (bottom, Fig- 
ure |5J). 
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FIG. 3: 90 and 99% allowed regions for the ID scan of 
LSND, overlayed on the 2D published allowed regions. Tri- 
angle points contain the 99% CL region, circle points contain 
the 90% CL region. 
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FIG. 4: Summed 2D Ax 2 compatibility grid from LSND, 
KARMEN2, MiniBooNE, Bugey. The star indicates the point 
of maximal compatibility (3.94%). 



LSND, MiniBooNE 

Table U presents results from the combination of LSND 
and MiniBooNE, not including the KARMEN2 result. 
The compatibility from the 2D analysis is actually lower 
than that found from the combination of LSND, KAR- 
MEN2, and MiniBooNE. KARMEN2 and MiniBooNE 
are complementary results; KARMEN2 has the most 
power in high Am 2 regions, while MiniBooNE is most 
sensitive to the lower Am 2 areas. The maximum com- 
patibility of the (LSND, MiniBooNE), and (LSND, KAR- 
MEN2, MiniBooNE) analyses is found in the low Am 2 
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FIG. 5: Maximum compatibility as a function of Am 2 for the 
ID LSND, KARMEN2, MiniBooNE, Bugey analysis. 



region where MiniBooNE has the most power. The in- 
clusion of KARMEN2 data adds two degrees of free- 
dom, but very little resolving power in this area of phase 
space. Figure [TO] illustrates this effect by overlaying the 
2D LSND 90 and 99% allowed regions with the Mini- 
BooNE and KARMEN2 2D 90% exclusion curves. 

MiniBooNE previously reported a 2% compatibility for 
the combination of LSND and MiniBooNE, found using 
the ID raster scan method [? ]. The prior result was 
calculated over a restricted Am 2 range (0.2 to 0.7 eV 2 ). 
The current analysis, when restricted to the same Am 2 
range, agrees with the previously published result. 



LSND, KARMEN2 

The maximum compatibility of LSND and KAR- 
MEN2, found using the 2D method, is 32.21%. This dif- 
fers from a previously reported compatibility of 64% [? 
]. There are two differences between the current anal- 
ysis and the study by Church et al.: the input LSND 
data set (this analysis utilizes the LSND decay-in-flight 
and decay-at-rest results, while the previous analysis only 
used the LSND decay-at-rest data) , and the method used 
to define and calculate the compatibility. Both analy- 
ses find a high compatibility between LSND and KAR- 
MEN2. 



CONCLUSIONS 

We present results on the compatibility of different 
combinations of four experiments which have searched 
for neutrino oscillations at the high Am 2 scale (> 0.0488 
eV 2 ). The LSND experiment has observed a significant 
excess of events; the other three experiments report 
null results and set limits on the oscillation parameter 
space. The compatibility has been calculated using 
both a 2D and a ID scan technique with the method 
of Reference [? ]. The remaining allowed regions 




2-D 90, 99% CL Exclusion Bands (KAR2, MB) 
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FIG. 6: Top: Summed Ax 2 compatibility of KARMEN2 and 
MiniBooNE using the 2D analysis. The star indicates the 
point of maximal compatibility (73.44%). The compatibility 
is limited by the boundaries of the analysis and may increase 
with a loosening of the grid range. Bottom: Exclusion bands 
(90%, 99%) found for the 2D KARMEN2 and MiniBooNE 
joint analysis. Values to the right of the lines are excluded 
at the 90, 99% CL. Triangle points form the 99% CL band, 
circle points form the 90% CL band. 



have been found for combinations resulting in greater 
than 10% compatibility. Results from the 2D scan 
indicate that LSND, KARMEN2, and MiniBooNE are 
25.36% compatible with having come from two-neutrino 
oscillations. However, the best fit point for this analysis 
is found in a region excluded by Bugey. (This point 
is also excluded by other reactor experiments such as 
Goesgen [? ], and Krasnoyarsk [? ].) The 2D scan from 
all four experiments including Bugey, in a limited sin 2 29 
region common to all experiments, finds they are only 
3.94% compatible with two-neutrino oscillations. This 
analysis does not take into consideration the absolute 
goodness of fit of each individual experiment at its own 



8 



Maximum 1-D Compatibility (KAR2, MB) 




best fit point, or any additional non-standard model 
effects such as CP violation or sterile neutrinos. 
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1-D 90, 99% CL Exclusion Curves (KAR2, MB) 
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FIG. 7: Top: Maximum compatibility as a function of Am 2 
for the ID KARMEN2 and MiniBooNE analysis. Bottom: 
Exclusion bands (90%, 99%) found for the ID KARMEN2, 
MiniBooNE joint analysis. Values to the right of the lines are 
excluded at the 90, 99% CL. Triangle points form the 99% CL 
band, circle points form the 90% CL band. These curves are 
very similar to those found using the 2D method in Figure [S] 
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FIG. 8: Top: Summed 2D Ax 2 compatibility grid from KAR- 
MEN2, MiniBooNE, and Bugey. The star indicates the point 
of maximal compatibility (27.37%). Bottom: Exclusion bands 
(90%, 99%) found for the 2D KARMEN2, MiniBooNE, and 
Bugey joint analysis. Values to the right of the lines are 
excluded at the 90, 99% CL. Triangle points form the 99% 
CL band, circle points form the 90% CL band. The vertical 
straight edge on the left indicates the lower sin 2 26 bound of 
0.01. 
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FIG. 9: Top: Maximum compatibility as a function of Am 2 
for the ID KARMEN2, MiniBooNE, and Bugey analysis. 
Bottom: Exclusion bands (90%, 99%) found for the ID KAR- 
MEN2, MiniBooNE, Bugey joint analysis. Values to the 
right of the lines are excluded at the 90, 99% CL. Triangle 
points form the 99% CL band, circle points form the 90% 
CL band. The vertical straight edge on the left indicates the 
lower sin 2 28 bound of 0.01. 
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MB, KARMEN2 90% CL Exclusion Bands 




FIG. 10: LSND 90, 99% CL allowed regions, overlayed with 
the KARMEN2 (triangles) and MiniBooNE (circles) 2D 90% 
CL exclusion bands. KARMEN2 is more powerful at exclud- 
ing high Am 2 values while MiniBooNE is more powerful in 
the low Am 2 region. 



